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ABSTRACT 

We have obtained HST NICMOS snapshot images of ten BALQSOs and four radio-quiet non-BALQSOs to 
search for extended narrow-line regions. For each object we obtained an image in the F160W filter and a narrow- 
band image in either Ha, H/3, or [O II]. We find no detections of extended narrow-line emission in the sample. We 
generate simulated images of QSOs with extended narrow-line regions to place limits on the amount of flux that 
may be present, and we conclude that BALQSOs cannot be more than several times brighter in extended emission 
than typical type 2 AGN. Spatially unresolved Ha emission is detected in all nine objects observed in this line, 
which we attribute to broad Ha emission. The five detected BALQSOs have derived Ha equivalent widths that 
are larger than those of the four non-BALQSOs by a factor of ~2 on average. This may point to an intrinsic 
difference between the rest-frame optical spectral properties of BALQSOs and non-BALQSOs. 

Subject headings: (galaxies:) quasars: absorption fines — (galaxies:) quasars: emission lines — (galaxies:) 
quasars: general 



1. INTRODUCTION 

Broad absorption line QSOs (BALQSOs) exfiibit broad ab- 
sorption troughs blue-shifted from the peaks of the correspond- 
ing broad emission lines, typically including Lya A1216, C IV 
A1549, Si IV A1397, N V A1240, and O VI A1034. They rep- 
resent a small fraction of optically-selected QSOs, about 9%, 
though the inferred fraction increases to about 12% when the 
luminosities are corrected for flux lost in the absorption troughs 
(Foltz et al. 1990). Additional effects such as scattering attenu- 
ation of the continuum (Goodrich 1997) or a tendency for BAL 
lines of sigfit to lie at large inclination angles with respect to the 
accretion disk (Krolik & Voit 1998) would make BALQSOs 
even more underrepresented in flux-limited surveys, and thus 
the true population fraction of BALQSOs could be significantly 
greater than 12%, perhaps as large as 50%. One of the most ba- 
sic and important problems concerning BALQSOs is the popu- 
lation fraction and how it relates to the location and geometry of 
the absorbing material. At one extreme, the BALQSOs could 
be a separate class of QSOs with a large, nearly unity, cov- 
ering fraction of absorbing clouds. Alternatively, in a unified 
view, the population fraction of BALQSOs is related to the typ- 
ical covering fraction of the absorbing gas. Through a detailed 
study of emission and absorption line profiles of BALQSOs, 
Hamann, Korista, & Morris (1993) conclude that the covering 
fraction of absorbing material is generally < 0.2. However, 
the corrections to the observed population fraction suggested 
by Goodrich (1997) and Ki-olik & Voit (1998) would imply 
larger covering fractions. If the typical covering fraction is in- 
deed much less than unity, then it seems plausible that BALQ- 
SOs and radio-quiet non-BALQSOs belong to the same popula- 
tion and differ only in their orientations relative to the observer, 
analogous to unified models of Seyfert 1 and 2 galaxies. Thus 
understanding the location and geometry of the absorbing gas 
in a BALQSO could be an important step in developing a full 



picture of the structure of all QSOs. 

Weymann et al. (1991) favor a geometry in which clouds ab- 
lated from a torus surrounding the central source are accelerated 
into an outflowing wind. Lines of sight lying near the plane of 
the torus for a given object would result in its classification as a 
BALQSO. One could also envision a model in which the high- 
velocity outflow was directed along the axis of the torus so that 
BALQSOs are viewed instead pole-on. Spectropolarimetric ob- 
servations with the Keck 10-m telescope (Goodrich & Miller 
1995; Cohen et al. 1995), however, favor a geometry in which 
the absorbers are in a flattened configuration around the nu- 
clear region associated with, or just outside, the broad emission- 
line region (BELR), consistent with the Weymann et al. (1991) 
model. The suggested geometry for BALQSOs strongly resem- 
bles the obscuration, reflection, orientation model that provides 
a unified view of Seyfert 1 and 2 galaxies (Antonucci 1993), 
although recent evidence suggests that some Seyfert 2 galax- 
ies are intrinsically distinct from Seyfert Is (Malkan, Gorjian, 
& Tam 1998). In the unified model, Seyfert Is are observed 
pole-on with a direct view of the central regions; the line of 
sight to Seyfert 2s intercepts the obscuring torus, but electron 
scattering provides a view of the nuclear regions. A corollary 
of this geometrical model is that shadowing by the torus can 
also collimate the ionizing radiation into cones that illuminate 
gas in the surrounding galaxy. In Seyfert 2s the illuminated gas 
shows extended, bipolar structures projected onto the sky, while 
in Seyfert Is the structures are expected to be more symmetric 
and compact (Mulchaey, Wilson, & Tsvetanov 1996; Schmitt & 
Kinney 1996). If this geometrical picture of BALQSOs is cor- 
rect, then BALQSOs, similar to Seyfert 2s, are preferentially 
viewed near the plane of the obscuring torus. If QSOs pos- 
sess extended narrow-line regions (NLRs) collimated by the 
torus as seen in Seyfert 2s, these extended regions will pref- 
erentially have larger projections on the plane of the sky than 
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non-BALQSOs. 

A significant difference between BALQSOs and Seyfert 2s 
is that in BALQSOs our line of sight to the nucleus and BELR 
is relatively clear (Weymann et al. 1991). A completely opaque 
torus along the line of sight, however, is not necessary to colli- 
mate the radiation. All that is required is material that is opaque 
to ionizing radiation. Some Seyfert Is, notably NGC 4151 and 
NGC 3516, show strong, broad UV absorption Unes (Penston 
1981; Uhich & Boisson 1983), opaque Lyman limits (Kriss et 
al. 1992, 1996), and extended, bipolar NLRs (Evans etal. 1993; 
Miyaji, Wilson, & Perez-Fournon 1992). Evans et al. (1993) 
and Kriss et al. (1995, 1996) suggest that the optically thick 
neutral hydrogen in these galaxies collimates the ionizing radi- 
ation. The same principal can work in BALQSOs. 

Extended narrow-line emission from QSOs has been read- 
ily observed and studied at lower redshift (z < 1), but has been 
mostly confined to radio-loud QSOs. The presence of extended 
[O III] emission has been found to be preferentially associated 
with steep-spectrum radio sources (Boroson & Oke 1984; Boro- 
son, Persson, & Oke 1985; Stockton & MacKenty 1987). Dur- 
ret et al. (1994) used integral field spectroscopy to study the ex- 
tended emission in greater detail and found that the morphology 
and velocity structure of the ionized gas was quite irregular and 
complex. Hes, Barthel, & Fosbury (1996) detected extended 
[O II] A3727 emission around several radio-loud QSOs. 

Following the analogy to Seyfert 2 galaxies, we undertook a 
snapshot program with HST NICMOS to search for extended 
Une emission that might be associated with an extended NLR 
in BALQSOs. No extended emission was detected in our pro- 
gram. In §2 and §3 we discuss our observations and analysis, 
and in §4 we discuss the impUcations of the lack of extended 
emission in our sample. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained images of 14 QSOs from 1998 March to 
November as an HST NICMOS snapshot program. Our full 
sample of 28 QSOs was chosen from Weymann et al. (1991) 
and includes 21 BALQSOs and 7 radio-quiet non-BALQSOs. 
Objects were selected by redshift such that the narrow emission 
hue of [O III] A5007 or A4959 (none observed), [O II] AA3726, 
3729 (4 objects observed). Ha A6563 (9 objects observed), or 
H/3 A4861 (1 object observed) was shifted to a wavelength 
within the bandpass of one of the narrow-band filter / camera 
combinations available on NICMOS. All data were acquired in 
MULTIACCUM mode. Each object was imaged in the narrow 
band for 703.94 s using the Step64 exposure time sequence and 
in the F160W filter with the same camera for 71.93 s using the 
Step8 exposure time sequence. The observation date and fil- 
ters used for each object are listed in Table 1. Unfortunately, 
the original redshift assigned to the BALQSO Q 0021-0213, 
Z = 2.293, which would have placed the [O III] emission line 
in the narrow band pass, was incorrect. The revised redshift, 
z = 2.348 (Hewett, Foltz, & Chaffee 1995), places the line out- 
side the band pass of F164N, and we do not consider this object 
as part of our observed sample. 

All of the images were processed with IRAF using the STS- 
DAS calibration pipeline task calnica along with the biaseq and 
pedsky tasks to remove the non-linear readout-to-readout bias 
levels and background pedestal, respectively. Some residual 
shading was evident in the NIC 2 images, which was removed 
with the IRAF task background, fitting the background perpen- 
dicular to the fast read direction. 



We identified a potential problem with the CRIDCALC step 
in the calnica pipeline that results in underestimates of the er- 
rors for MULTIACCUM data. To produce the final science and 
error images, calnica fits a line to the counts vs. time data pro- 
vided by the many MULTIACCUM images and derives the fi- 
nal count rate and error per pixel based on this fit. For data 
that are read-noise dominated, this is a vaUd approximation. 
However, the algorithm does not take into account the fact that 
the Poissonian errors from read to read in the signal from any 
sources are correlated, and thus underestimates the errors in re- 
gions of high signal. We generated our own error images using 
the following formula: 



(»Mgc|*EXPTIME*ADCGAIN , ■nt/^tctttt T32 
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error=— , (1) 

EXPTIME*ADCGAIN 

where image is the final science image produced by calnica in 
units of DN sec-^ FLATFILE and NOISFILE are the particular 
flatfield image and RMS read noise image in electrons used by 
the pipeline, EXPTIME is the exposure time in seconds, and 
ADCGAIN is the analog-to-digital conversion gain in e~ DIST^ . 
This method overestimates the read noise contribution, since 
this will in general be smaller than the noise per read repre- 
sented by the NOISFILE due to the multiple readouts used in 
MULTIACCUM mode. We do not include the error from the 
dark current contribution. In order to check that ignoring the 
dark current in this calculation is justifiable, we estimate the 
typical dark current contribution to our images using typical 
values for the dark current contribution as described in Skinner 
& Bergeron (1997). The dark current consists of three basic 
components: the true linear dark current, amplifier glow from 
infrared radiation emitted by the readout amplifiers, and shad- 
ing due to time dependence of the bias level during readout. The 
shading component is noiseless. The Unear dark current is only 
^ 0.05 s"' pixel"', so even with a 700 s exposure, the error 
contribution from the linear dark current is only ~ 6 e~ pixel"' 
and is thus negligible compared to the read noise, which is typ- 
ically ~ 30— 40 e~ pixel"'. The largest error contribution from 
the dark current is from amplifier glow, which can contribute 
2-3 DN pixel"' per read to the signal. The gain for our data is 
5.0 e" DN"' for cameras 1 and 2, 6.5 e" DN"' for camera 3. 
Given that there are 20 readouts in the Step64 exposure time 
sequence, and assuming 3 DN pixel"' per read from amplifier 
glow and a gain of 5.0, the amplifier glow would contribute 
20 X 3 X 5 = 300 pixel"' to the signal in the final readout, 
for an error contribution of ~ 17 e" pixel"'. The contribution is 
thus not insignificant but still small compared to the read noise 
that we include, and we ignore it for simplicity. Overall we 
have found that the error images produced by this method yield 
a good approximation of the error in regions of low signal and, 
by design, a much better representation of the Poissonian er- 
rors. We thus consider it to be a more reliable indicator of the 
true uncertainty for our particular data sets and use it in our 
analysis. 

3. ANALYSIS 

We use two different methods to determine the presence and 
magnitude of any emission-line flux. In the first method, we use 
the broad-band image to determine the continuum contribution 
to the count rate in the narrow-band image in order to calcu- 
late the net emission-Une flux. This method ignores the spatial 
distribution of the line emission. In the second method, we 
measure the residual count rate in PSF- subtracted images. To 
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quantify limits on the amount of extended line emission from 
these residuals, we create models of an extended NLR com- 
bined with a point source as well as PSF- subtracted models to 
compare with the real data. 

In both of these methods we utilize model PSFs created with 
the HSTPSF generation program Tiny Tim (Krist 1993) version 
5.0. All the PSFs are calculated for the observed position on 
the camera array and for the appropriate observation date. The 
PSFs are subsampled by a factor of 10 for better accuracy in 
sub-pixel shifting, used for the PSF subtraction. For simplicity, 
we use monochromatic PSFs, using for the representative wave- 
length the effective wavelength of the bandpass for the narrow- 
band PSFs and 1600 nm for the broad-band PSFs. Although 
in general one might not expect a monochromatic PSF to be a 
very good match for F160W data, we use the broad-band PSFs 
only to obtain aperture correction factors, and our tests show 
that for our apertures these factors vary by < 1 % throughout 
the bandpass (~ 1400- 1800 nm), so they are suitable for this 
purpose. 

Using aperture photometry, we measure the brightness of 
each QSO in the broad-band and narrow-band images. We 
choose aperture radii of an integral number of pixels such that 
the aperture encloses at least 90% of the light based on model 
PSFs. The resulting apertures enclose the Airy disc and the first 
bright Airy ring and have radii between 0"5 and 0"8 in angular 
extent, depending on the camera and filter being used. This size 
aperture is suitable since it is large enough to be insensitive to 
minor variations in the PSF and to enclose a significant amount 
of the extended emission, if it is present, but not so large as to 
significantly increase the uncertainty of the measurement due 
to additional background and read noise. These count rates are 
then multiplied by an appropriate aperture correction factor as 
calculated from the model PSFs. Table 2 lists the raw and cor- 
rected count rates with errors as obtained by this method. The 
errors attributed to the measurements are simply the quadrature 
sums of the errors of each pixel included in the aperture, calcu- 
lated as described in §2. We also list the F160W magnitudes of 
each object in AB mags. 

To derive estimates of the emission-line flux from these data, 
we use the STSDAS package synphot to calculate the expected 
ratio of broad-band to narrow-band continuum count rates for 
each object assuming a power-law continuum oc i/-0-3±0-6 
The power-law index and RMS deviation were estimated by in- 
spection from the statistical data of Francis et al. (1991), which 
covers the continuum wavelength range of interest. This allows 
us to scale the broad-band fluxes to remove the continuum con- 
tribution in the narrow band. We use synphot to shift a model 
line proflle to the proper wavelength and calculate conversion 
factors from net emission-line count rates to flux for each ob- 
ject. For [O IIJ and [O III] we model the line as a Gaussian 
with a FWHM of 400 km s"^ Ha and H/3 are both modeled 
by the broad H/3 profile extracted from the Francis et al. (1991) 
composite. (Ha is outside the wavelength coverage of the com- 
posite.) As applied to Ha, this model profile does not include 
any contribution from the narrow lines [N II] AA6548, 6583, but 
in broad-line objects these lines make a very small contribution 
to the emission-line flux compared to broad Ha, as can be seen 
for example in Jackson & Browne (1991). The resulting con- 
version factors are shown in column 5 of Table 3. 

The model Balmer profile has a FWHM of 2200 km s~' and 
a full width at zero emission of 7700 km s"\ whereas the nar- 
row bandpasses are generally rectangular in shape with a width 



of 2500-3000 km s~'. The conversion factors thus represent 
a conversion from count rates to total integrated emission-line 
flux, including for Ha and H/3 portions of the broad component 
that lie outside the range of the narrow-band filter. The con- 
versions are therefore sensitive to the assumed profile shape. 
The model fine profile FWHM of 2200 km s~' is fairly narrow 
for a broad line. The seventeen radio-quiet QSOs in the z ^ 2 
Mcintosh et al. (1999) sample have a mean H/3 FWHM of 5 100 
km s~\ while the z < 0.5 Boroson & Green (1992) sample has 
a mean H ^ FWHM of 3800 km s"'. The BALQSOs and non- 
BALQSOs in the Weymann et al. (1991) study have mean C IV 
half- widths at half maximum of 2200 and 2500 km s~\ respec- 
tively. To illustrate the effect the line profile has on our results, 
we recalculate the conversion factors for the broad lines assum- 
ing a Gaussian profile with a FWHM of 5000 km s~' . These 
conversion factors are fisted in column 6 of Table 3, and the 
ratio of these values to those derived using the Francis et al. 
(1991) profile are listed in column 7. As the table shows, the 
conversion factors increase by various amounts up to ~ 40%, 
which would result in corresponding increases in the fluxes and 
equivalent widths we derive. The conversions are also sensi- 
tive to the position in wavelength space of the line profile in 
the narrow-band filter, since for a given line flux the count rate 
will be higher if the fine peak is near the center of the band- 
pass than if it is near the edge. To iUustrate this, we calculate 
the conversion factor that would be used if the peak of the fine 
were exactly at the effective wavelength of the narrow-band fil- 
ter, using the Francis et al. (1991) profile. These values are 
fisted in column 8 of Table 3. The ratio of the conversion fac- 
tors using the true wavelengths to those assuming the emission 
lines are centered are listed in column 9. It is perhaps useful 
to consider column 9 as a correction factor for the fact that the 
fine is not centered in the bandpass. As one can see by com- 
paring these values to the emission-line peak wavelengths and 
filter effective wavelengths in columns 3 and 4, we are correct- 
ing significantly more for flux outside the bandpass when the 
fine is not centered, a point which is relevant to the comparison 
of the BALQSOs and non-BALQSOs in our sample discussed 
in §4. 

The emission-line fluxes, calculated using the conversion 
factors from column 5 of Table 3, are shown in Table 4. The 
errors include an assumption of an uncorrelated 2% uncer- 
tainty in the absolute throughput of each of the filters, although 
this makes a negUgible contribution in comparison to the un- 
certainty in the continuum slope. (Information on NICMOS 
photometric calibrations can be found at the NICMOS web 
site, http://www.stsci.edu/cgi-bin/nicmos.) The errors do not 
include any uncertainty in the profile shape. We also Ust in 
the table the corresponding equivalent widths of the emission 
lines in the observer and object rest frames. It is interesting 
that of the nine objects we observed in Ha, the five BALQSOs 
have systematically larger equivalent widths than the four non- 
BALQSOs. We discuss fliis further in §4. 

Using aperture photometry to measure the emission-line flux 
is not particularly sensitive as it combines several sources of un- 
certainty, most notably the throughputs of the filters, the QSO 
continuum spectral shape, and the profiles of the emission lines. 
As a result, the only solid detections of any emission-line flux 
by this method are of Ha and H/3, most likely dominated by 
the broad component. No forbidden fines typically seen in ex- 
tended NLRs, such as [O n] AA3726, 3729, are detected. 

Although the observations permit us to measure the strength 
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of the emission lines by this method, the primary purpose of 
the observations was to search for extended emission. None 
of the objects show any obvious extended emission from a 
simple visual inspection of the images (see Figures 1-3). To 
probe more sensitively for extended narrow-line emission, we 
need to remove the point-source contribution. We experimented 
with subtracting scaled broad-band images from the narrow- 
band images to create a continuum-subtracted image as is typ- 
ically done for extended sources, but because the emission is 
(at least mostly) point-like, the resulting image is dominated 
by the residuals from the variation in PSFs between bands. In- 
stead, we create PSF- subtracted narrow-band images, using the 
model PSFs generated individually for each object with Tiny 
Tim as discussed previously. 

We use a simple iterative process to find the best-fit subtrac- 
tion by minimizing the RMS deviation of the fit residuals. The 
position and normalization of the point source are varied as free 
parameters in the fitting process. The PSF subtractions worked 
quite well in general, particularly for cameras 1 and 2, although 
the camera 3 subtractions did leave some residuals. The PSF 
subtracted images are shown in the second row of Figure 1 for 
the NIC 1 images. Figure 2 for NIC 2, and Figure 3 for NIC 
3. We measure the residual count rate after PSF subtraction in 
the same fixed aperture as before. These count rates are con- 
verted to flux, assuming the residual counts are due to a 400 
kms"' FWHM emission Une at the expected wavelength for 
all lines. The results are shown in column 2 of Table 5. The 
quoted errors only include the individual pixel noise, so they do 
not include any uncertainty in the fits or the PSFs. Since we 
find residual fluxes both positive and negative well above this 
noise level, there are certainly additional sources of uncertainty. 
From inspection of the images, we conclude that all of the resid- 
ual fluxes can be explained by either imperfect PSF subtractions 
or by variations in the background, probably associated with 
flatfield features not entirely removed by the pedestal correc- 
tion. 

As a final check on whether any extended emission is present 
in our images and to set upper Umits on its absence, we create 
models to simulate the expected appearance of the extended 
emission in the image and the residuals from the PSF sub- 
traction. To model Ihc shape of the extended emission we 
use an appropriately scaled continuum- subtracted WFPC2 im- 
age in [O III] of NGC 1068 (Dressel et al. 1998). We as- 
sume the same size scaling for the extended emission of each 
model. The angular scale at the distance of NGC 1068 is 72 
pc arcsec"' (Bland-Hawthorn et al. 1997). Assuming fio = 1 
and flo = 100 km s"' Mpc"' , the angular scale at z = 2, approx- 
imately the mean redshift of our sample, is 4. 1 kpc arcsec"' . 
For h = 0.75, the angular scale is 5.5 kpc arcsec"^ To scale 
NGC 1068 up to the luminosity of a QSO, we assume that the 
NLR is optically thin and radiation bounded. Thus the physical 
size should scale roughly as the square root of the intensity of 
the ionizing radiation, which should scale roughly with bolo- 
metric luminosity. The bolometric luminosity of NGC 1068 
is Lhoi^Nccw6i ^ 4.2 X lO'*'* erg s"' (Pier et al. 1994), whereas 
a typical QSO in our sample has a bolometric luminosity of 
Lboi,QSO ^ 2 X 10"** ergs"'. We thus expect the NLR to be 
~ ^2 X 10''*'/4.2 X 10'*^ = 6.9 times larger in hnear physical 
extent. Including the relative angular scales due to distance as 
described above, the overall size scaling to be appUed to the 
NGC 1068 image to make it resemble a hypothetical BALQSO 
is 0.09. We scale the NGC 1068 image by this amount and re- 



sample it to match the pixel size of the various NICMOS cam- 
eras. We normalize the flux for each model such that the total 
flux of the extended emission is some particular fraction of the 
aperture-corrected flux. We then convolve the model with the 
appropriate PSF and place it onto an empty section of the orig- 
inal image. A point source is added to the model, representing 
the remainder of the flux, such that the total flux of the point 
source plus extended component always equals the aperture- 
corrected value. This is then fit and subtracted using the same 
method as for the real data. We repeat this process for various 
fractions of extended emission relative to total flux in order to 
find at what minimum fraction the extended emission becomes 
detectable. We judge the extended emission to be detectable 
if some structure of several contiguous pixels above (or below) 
the background noise is visually apparent in the residuals. This 
process allows us to place an upper limit on the observed ex- 
tended emission, albeit a subjective one. 

The models with the minimum detectable fraction of ex- 
tended emission are shown in row 3 of Figures 1-3. The PSF- 
subtracted models are shown in row 4 of the figures. Columns 

3 and 4 of Table 5 list the narrow emission-line fluxes corre- 
sponding to the minimum detectable flux models and the resid- 
uals from the PSF subtractions of these models. These residuals 
were measured in the same fixed aperture as the other flux mea- 
surements and can thus be directly compared to the values from 
the real data in column 2. The residuals from the models are 
in most cases clearly larger than what was measured from the 
real data, consistent with our conclusion that we do not detect 
any extended flux. It is evident from comparing columns 3 and 

4 of Table 5 that the fit requirement of minimizing the RMS 
deviation of the residuals results in a substantial amount of the 
extended flux being removed by the PSF subtraction, as much 
as - 90%. 

Our tests using the models suggest that the presence of ex- 
tended emission would result in excess flux in the first dark Airy 
ring as well as an over-subtracted core in the PSF-subtracted 
images, as can be seen in row 4 of Figures 1-3. The camera 1 
and 2 data clearly show no such features in the PSF-subtracted 
images. For the only object observed with camera 1 or 2 with 
any significant residuals, Q 2358H-0216, the residuals are along 
the first bright Airy ring and in the core. This is more indicative 
of an imperfect PSF subtraction than of extended emission. The 
camera 3 data are more difficult to judge. Although in princi- 
ple camera 3 is most sensitive for our purpose due to the larger 
pixel scale and correspondingly lower read noise per unit an- 
gular area, the undersampling of the PSF makes the fits less 
certain and makes it difficult to determine if there is excess flux 
in the first dark Airy ring, which has a width of only ~ 0.5 pix- 
els. Since the residuals from the camera 3 subtractions do not 
resemble closely those expected from the models, we conclude 
that the residuals are due to problems with the subtraction and 
are not indicative of extended emission. 

4. DISCUSSION 

To provide a context for the limits we have placed on the 
amount of extended emission, we use two separate methods 
to estimate the amount of extended emission we expected to 
observe. For the first estimate we assume that the equivalent 
widths of the narrow forbidden emission lines should be the 
same as for the composite spectrum of Francis et al. (1991). 
The equivalent widths for the composite are 15 A for [O III] 
A5007 and 1.9 A for [O II] A3727. To estimate narrow-line 
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equivalent widths for the Balmer lines, we assume flux ratios 
of 10 for [O III] A5007 / narrow H/3 and 3 for narrow Ha / 
narrow H/3. These ratios are typical for Seyfert 2s and narrow- 
line radio galaxies (Koski 1978). Assuming a continuum shape 
/v oc i^"" -', we derive expected rest-frame equivalent widths of 
~7.5 A for narrow Ha and ~ 1 .5 A for narrow H/3. These equiv- 
alent widths are used to estimate narrow-line fluxes for each 
object, listed in column 5 of Table 5. 

For the second estimate of the expected narrow-line flux we 
scale the observed flux of the emission lines in NGC 1068 
(Koski 1978) to the luminosity and distance of the QSOs in 
our sample. We assume that the flux of the lines scales linearly 
with bolometric luminosity. The bolometric luminosities of the 
QSOs are estimated as XF\ as observed in the F160W filter. 
The luminosity distance for each QSO is calculated from the 
redshift assunoing = 1 and h = 0.75. The Une fluxes scaled 
from NGC 1068 are listed in column 7 of Table 5. 

Columns 6 and 8 of Table 5 list the ratios of these fluxes to 
the minimum detectable flux in column 3 determined from the 
models. A ratio greater than unity would indicate that we would 
expect to have detected extended emission. However, only the 
NGC 1068 scaling for PHL 5200 exceeds unity. This object 
had the poorest PSF subtraction, so it is difficult to ascertain 
the existence of extended emission. Depending on the scaUng, 
roughly half of the objects are above the 10% level, while some 
of the others are much smaller, well below our detection thresh- 
old. If the extended emission is at the level expected by these 
scaling arguments, it is not particularly surprising that we did 
not detect it. However, if the emission were only a few times 
stronger than what we expect based on these simple scaling 
arguments, then we should have been able to detect extended 
emission in some cases. 

It is unfortunate that none of the six objects on the target list 
for our snapshot program matched for [O III] were observed, 
since the [O III] lines are expected to be the strongest emis- 
sion lines in the NLR. For our scaling argument above, we as- 
sumed an [O in] A5007 equivalent width of 15 A. Mcintosh 
et al. (1999) use infrared spectroscopy to study the rest-frame 
optical emission-line properties of a sample of QSOs with sim- 
ilar redshifts and luminosities to our sample. They find a mean 
equivalent width for [O III] A5007 of 8.1 A for seven BALQ- 
SOs and 11.6 A for ten radio-quiet non-BALQSOs, but given 
the uncertainty and the small sample size the difference is not 
significant. These equivalent widths are somewhat lower than 
the mean of 23 A for the 87 z < 0.5 PG QSOs observed by Boro- 
son & Green (1992). The low-redshift IR-selected QSO sam- 
ple of Boroson & Meyers (1992) has even higher mean [O III] 
equivalent width of 46 A for the 14 objects which were not low- 
ionization BALQSOs. Given these results, with even a conser- 
vative estimate of the strength of [O III], we would expect to 
observe [O III] line fluxes a factor of 3-5 greater than what we 
list in Table 5. We therefore might reasonably have expected to 
see extended emission in [O III] in 1-2 objects had they been 
observed. 

The observed difference in Ha equivalent widths be- 
tween BALQSOs and non-BALQSOs merits further discus- 
sion. The mean rest-frame Ha equivalent width for the BALQ- 
SOs is 390 ± 30 A, while the mean equivalent width for 
the non-BALQSOs is 200 ± 20 A, corresponding to a ratio 
(EWh„(BAL)) / (EWHa(non - BAL)) = 1 .9 ± 0.2. Even more 
striking than the ratio of the means is that the difference is 



systematic — all of the BALQSOs have larger Ha equivalent 
widths than all of the non-BALQSOs. A two-sample K-S test 
yields that the BALQSO and non-BALQSO data were drawn 
from separate populations with 99.3% confidence, although 
given that the sample groups are small and certainly incom- 
plete, one should approach this number with caution. However, 
the data are suggestive of a possible real systematic difference 
between the two groups. 

We have considered the possibility that the difference in 
equivalent widths could be due to some systematic difference 
in the observations. The observations were performed over 
roughly the same time period and with the same filters and cam- 
eras for both groups. However, there is a systematic difference 
in the location of the Ha line in the filters. For all of the non- 
BALQSOs, the expected location of the peak of the broad Ha 
emission Une lies near the center of the filter bandpass, while 
for the BALQSOs the peak of the line lies near the edge of 
the bandpass, resulting in larger corrections for flux outside the 
bandpass. This can be seen by comparing the values in col- 
umn 9 of Table 3 for the BALQSOs and non-BALQSOs. For 
three of the BALQSOs the emission line center is longward 
of the center of the filter, while for the other two it is short- 
ward. The systematic difference in the location of the lines in 
the filters renders the effect of making the relative Ha equiv- 
alent widths of the BALQSOs and non-BALQSOs sensitive to 
the assumed profile shape. We find that if we assume a broader 
line profile, the disparity between the two groups diminishes 
but does not disappear. If we use the conversion factors from 
column 6 of Table 3, calculated assuming a Gaussian profile 
with a FWHM of 5000 km s"\ the ratio of equivalent widths 
becomes (EWHa(BAL))/(EWHa(non-BAL)) = 1.6±0.2. 

There are several possible systematic differences between 
BALQSOs and non-BALQSOs that could explain our observa- 
tions. The most straightforward explanation is that there is sim- 
ply more Ha flux in BALQSOs than non-BALQSOs. Another 
possibility is that the Ha lines could have the same flux on aver- 
age but be systematically much broader in the non-BALQSOs. 
Either one of these would be extremely surprising, however, 
since broad emission lines in the UV exhibit no large system- 
atic differences (Weymann et al. 1991). Also, Mcintosh et al. 
(1999) do not find large differences in H/3 emission between 
BALQSOs and non-BALQSOs. An even more intriguing pos- 
sibility is that the continua are different in BALQSOs and non- 
BALQSOs, since in all cases we have estimated the continuum 
in the narrow-band filter from shorter wavelength F160W im- 
ages. If the continuum shape were redder than we have as- 
sumed in BALQSOs, we would underestimate the continuum 
in the narrow-band images of the BALQSOs, resulting in larger 
fluxes and equivalent widths. The implied continuum shape 
in the rest-frame optical for BALQSOs that would account for 
the systematically larger Ha equivalent widths has an effective 
power-law index a^ = -1.9, compared to the typical a^ = -0.3 
that we assumed. This amount of reddening is comparable 
with that observed in low-ionization BALQSOs (Weymann et 
al. 1991) and attributed to intrinsic dust extinction (Sprayberry 
& Foltz 1992). In fact, two of the BALQSOs in our sample, 
PHL 5200 and Q 2358-1-0216, are low-ionization BALQSOs, so 
it is not surprising that we derive larger Ha equivalent widths 
for these objects. However, the fact that the other three BALQ- 
SOs follow the same trend is interesting, since there is gen- 
eraUy Uttle evidence from rest-frame UV spectra for intrinsic 
dust extinction in high-ionization BALQSOs (Weymann et al. 
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1991). To investigate whether there is an overall reddening of 
the continuum in our sample, we used our fluxes in the F160W 
band for each object along with the fluxes at ~ 2100 A in the 
rest frame calculated from optical spectra by Weymann et al. 
(1991) to derive (R-H) color indices. For the two QSOs in 
our sample not in Weymann et al. (1991), Q 0054+0200 and 
Q 2358+0216, we estimate the optical flux from the Bj magni- 
tudes in the LBQS catalog (Hewett, Foltz, & Chaffee 1995). Of 
the objects observed in Ha, we find that both the BALQSO and 
non-BALQSO groups have {(R-H)) = 1 .95 with a - 0.8 mag, 
suggesting that there is no significant systematic difference in 
the continuum properties in our sample in the rest-frame 2100- 
5500 A spectral region. Near-infrared (rest-frame optical) spec- 
troscopy of BALQSOs is necessary to resolve this issue. 

5. SUMMARY 

We have obtained HST NICMOS snapshot images of ten 
BALQSOs and four non-BALQSOs to search for extended 
narrow-line emission. Each object was imaged in both the 
F160W filter and a narrow-band filter corresponding to Ha, 
H/3, or [O II]. Using aperture photometry, we derive estimates 
of the emission-Une fluxes and equivalent widths. To look 
for extended emission, we use model PSFs to generate PSF- 
subtracted narrow-band images. We measure the residual flux 
in these images and compared them to PSF-subtracted images 
of simulated data. Our most important results are: 

1. We do not detect extended emission in any of the 

images. Some of the images have PSF-subtraction 
residuals, but based on comparisons with simulated 
data we conclude that the residuals are the result of 
imperfect PSF subtraction and are not evidence for 
extended emission. 

2. We find that for the nine objects we observed in Ha, 
we derive systematically larger equivalent widths for 
flie five BALQSOs than for the four non-BALQSOs, 
by a factor of 1.9 on average. Given the small sample 
size, it is not certain whether this is a real effect or 
just a statistical aberration. If it is real, it could be 
evidence either for an intrinsic difference between 
BALQSOs and non-BALQSOs in Ha emission, or 
for systematically redder rest-frame optical continua 
in BALQSOs, an effect which has been observed 

in low-ionization BALQSOs. If either of these 
differences were real, it would stand in stark contrast 
to studies of rest-frame UV spectra which find little 
difference between the continuum and emission-Une 
properties of high-ionization BALQSOs and radio-quiet 
non-BALQSOs. 
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Table 1 
Observations 



Broad-Band Narrow-Band 



Object 


Redshift 


BAL? 


Date 


Camera 


Filter 


Filter 


PHL 5200 


1.981 


Y 


1998 June 16 


3 


F160W 


F196N 


Ki 0021-0213 


'J AO 

2.348 


Y 


1998 October 10 


1 


F160W 


F164N 


Q 0054H-0200 


1.872 


Y 


1998 August 3 


2 


F160W 


F187N 


Q 0932H-5010 


1.914 


Y 


1998 September 13 


1 


F160W 


F108N 


Q 1029-0125 


2.029 


Y 


1998 June 20 


3 


F160W 


F200N 


Q 1146H-0207 


2.054 


N 


1998 March 15 


3 


F160W 


F200N 


Q 1225H-1512 


1.990 


N 


1998 June 14 


3 


F160W 


F196N 


Q 1333H-2840 


1.908 


Y 


1998 November 10 


1 


F160W 


F108N 


Q 1433-0025 


2.042 


N 


1998 May 8 


3 


F160W 


F200N 


Q 1439H-0047 


1.857 


N 


1998 July 10 


2 


F160W 


F187N 


Q 2154-2005 


2.035 


Y 


1998 October 26 


1 


F160W 


F113N 


Q 2358-1-0216 


1.872 


Y 


1998 August 16 


2 


F160W 


F187N 


UM 139 


2.029 


Y 


1998 October 18 


1 


F160W 


F113N 


UM 141 


2.909 


Y 


1998 September 27 


2 


F160W 


F190N 


UM253 


2.253 


Y 


1998 October 31 


2 


F160W 


F212N 



*Because of an incorrect redshift, no emission line is in the F164N bandpass, so we ignore this object in our analysis. 



Table 2 

APLRl'LiRh PlIO rOMh TR^' 



Broad-Band Count Rates* Narrow-Band Count Rates* 



Object 


Measured 


Corrected 


Measured 


Corrected 


ABfK^Qw 


PHL 5200 


249.7 ±0.8 


258.4±0.8 


62.9±0.1 


71.8±0.1 


16.821 ±0.003 


Q 0021-0213 


105.1 ±3.1 


115.8±3.4 


5.1 ±0.3 


5.6±0.4 


17.793 ±0.032 


Q 0054H-0200 


75.9±1.4 


84.4 ±1.5 


13.2±0.2 


14.6 ±0.2 


18.294 ±0.019 


Q 0932H-5010 


273.8 ±2.5 


301.6±2.7 


4.1±0.2 


4.5 ±0.2 


16.754 ±0.010 


Q 1029-0125 


113.0±0.7 


116.9±0.7 


22.1±0.1 


25.3±0.1 


17.682 ±0.006 


Q 1146H-0207 


123.3 ±0.7 


127.6±0.7 


14.8±0.1 


17.0±0.1 


17.505 ±0.006 


Q 1225H-1512 


74.4 ±1.2 


77.0 ±1.2 


16.4±0.1 


18.7±0.1 


18.136±0.017 


Q 1333H-2840 


96.7 ±3.5 


106.5 ±3.8 


1.5±0.3 


1.6 ±0.3 


17.884 ±0.039 


Q 1433-0025 


64.8 ±0.6 


67.1 ±0.6 


11.0±0.1 


12.1 ±0.1 


18.204 ±0.010 


Q 1439H-0047 


68.1±1.5 


75.7 ±1.7 


10.7±0.2 


12.1 ±0.2 


18.413 ±0.024 


Q 2154-2005 


155.2±3.2 


171.0±3.5 


3.9±0.3 


4.2±0.3 


17.370 ±0.022 


Q 2358H-0216 


214.7 ±1.4 


238.7 ±1.5 


35.2±0.2 


39.0 ±0.2 


17.165 ±0.007 


UM 139 


78.5±3.0 


86.5 ±3.3 


1.3±0.3 


1.4 ±0.3 


18.110±0.041 


UM 141 


207.4 ±1.4 


230.6 ±1.5 


13.8±0.1 


15.3 ±0.2 


17.203 ±0.007 


UM253 


97.4 ±1.4 


108.3 ±1.5 


21.6±0.2 


23.8 ±0.2 


18.023 ±0.015 



*DN s-i 
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Table 3 

Emission Line Flux Conversion Factors 





Emission 


Redshifted 


Filter 












Object 


Line 




A.//(A) 


A 


B 


B/A 


C 


A/C 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


PHL 5200 


Ha 


19569 


19639 


1.986 


2.430 


1.22 


1.017 


1.95 


Q 0054+0200 


Ha 


18854 


18740 


2.976 


3.259 


1.10 


1.805 


1.65 


Q 0932+5010 


[O II] 


10865 


10816 


13.40 






7.592 


1.77 


Q 1029-0125 


Ha 


19884 


19975 


2.147 


2.480 


1.10 


1.584 


1.36 


Q 1146+0207 


Ha 


20048 


19975 


1.823 


2.354 


1.29 


1.584 


1.15 


Q 1225+1512 


Ha 


19628 


19639 


1.621 


2.193 


1.35 


1.017 


1.59 


Q 1333+2840 


[Oil] 


10842 


10816 


7.665 






7.592 


1.01 


Q 1433-0025 


Ha 


19970 


19975 


1.589 


2.111 


1.33 


1.584 


1.00 


Q 1439+0047 


Ha 


18755 


18740 


1.807 


2.438 


1.35 


1.805 


1.00 


Q 2154-2005 


[Oil] 


11316 


11298 


5.880 






6.110 


0.96 


Q 2358+0216 


Ha 


18854 


18740 


2.976 


3.259 


1.10 


1.805 


1.65 


UM 139 


[Oil] 


11294 


11298 


6.279 






6.110 


1.03 


UM 141 


up 


19008 


19004 


1.760 


2.420 


1.38 


1.759 


1.00 


UM253 


Ha 


21355 


21213 


2.722 


2.757 


1.01 


1.391 


1.96 



Col. (3) — ^Expected wavelength in vacuum of the peak of the redshifted emission hne. 
Col. (4) — ^Effective wavelength of the narrow band filter. 

Col. (5) — Case A: Conversion factors from count rate to integrated emission-line flux assuming a 400 km s"' FWHM Gaussian for 
the narrow lines and the H/3 profile extracted from Francis et al. (1991) for the broad lines. Units are 10"' erg cm"^ s"' (DN s^')^'. 

Col. (6) — Case B: For the broad fines only, conversion factors assuming a 5000 km s~' FWHM Gaussian. Units are 10"'^ 
erg cm"^ s"' (DN s"')~^ 

Col. (7)— Ratio of B to A. 

Col. (8) — Case C: Conversion factors that would be used if the emission fines were centered in the narrow band filters, assunung 
the same profiles as A. Units are 10"^^ erg cm~'^ s~' (DN s~')~^ 

Col. (9)— Ratio of A to C. 
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Table 4 

Emission Line Data 



Narrow-Band Count Rates^ Equivalent Width (A) 

Scaled Net Integrated 

Object Line Continuum Emission Line Line Flux'' Observed Rest Frame 



PHL 5200 


Ha 


^^•^-2.0 


56.6!?:« 


1131 


1300!i« 


438!^^ 


Q 0054+0200 


Ha 




10.19!Oio 


30.3!J:^ 


1050!j^o 


364^6 


Q 0932+5010 


[Oil] 


A C!+0.9 




0+16 
^-13 


8+54 
°-43 


0+18 
•'-IS 


Q 1029-0125 


Ha 


'•^-1.0 




38.9!|i 


1030!?io 


TaQ+69 
■'•'^-60 


Q 1 146+0207 


Ha 


7.8^1? 


Q O+l.l 

^•^-1.0 


16.7^2-2 


400i''« 


133lt 


Q 1225+1512 


Ha 


4 C4+0.54 


14 14+0.62 


22.9+}-^ 


890!|» 


299^53 


Q 1333+2840 


[Oil] 




io+"-5" 

U.iU_„44 


0.8!3-f 


7+38 


0+13 
-'-11 


Q 1433-0025 


Ha 




7 QA+0.60 
' -^"-0.52 


12 7+'-0 


580+130 


191+42 


Q 1439+0047 


Ha 


■-'•^0-0.42 


7 87+0.45 


14.2i| 


ccn+90 


I9l!i^ 


Q 2154-2005 


[Oil] 


3.03t^« 


1 19+0.77 






14!1| 


Q 2358+0216 


Ha 


12.5!ji 


26.5!|i 


79^5 


960!li0 


335+52 


UM 139 


[Oil] 


^•■^■'-0.37 


-0 1 7+0-46 
^•^ '-0.40 


-1 1+2' 


_14+33 




UM 141 


H/3 




3 5+1-4 


6.1!|1 


79!^^ 


20+^1 


UM253 


Ha 




17.2^{i 


46.8^i:^ 


1550^^10 


480^120 



''DN s-i 

''10"'^ergcm-2 s"' 



Table 5 
PSF Subtraction Residuals 











Composite QSO Scaling 


NGC 1068 ScaUng 




Residual 


Minimum 


Model 


Expected 




Expected 




Object 


Flu.K^' 


Detectable Flux" 


Residual Flux" 


Line Flux" 


Ratio'' 


Line Flux" 


Ratio'" 


PHL 5200 


2.8 ±0.1 


4.9 


0.4 


1.9 


0.39 


6.0 


1.21 


Q 0054+0200 


-0.5 ±1.1 


33 


7.7 


0.6 


0.02 


1.5 


0.05 


Q 0932+5010 


-4.8 ±2.7 


30 


19.6 


1.7 


0.06 


LI 


0.04 


Q 1029-0125 


0.1±0.1 


6 


0.6 


0.9 


0.16 


2.7 


0.49 


Q 1146+0207 


-0.1 ±0.1 


10 


1.1 


0.9 


0.09 


3.2 


0.30 


Q 1225+1512 


0.5 ±0.1 


13 


1.3 


0.6 


0.05 


1.8 


0.14 


Q 1333+2840 


1.8 ±2.3 


16 


5.1 


0.6 


0.04 


0.4 


0.02 


Q 1433-0025 


0.5 ±0.1 


7 


0.6 


0.5 


0.07 


1.7 


0.23 


Q 1439+0047 


0.9 ±0.3 


5 


0.9 


0.6 


0.12 


1.4 


0.30 


Q 2154-2005 


4.5 ±1.8 


15 


6.2 


0.9 


0.06 


0.6 


0.04 


Q 2358+0216 


-2.2±1.1 


22 


6.4 


1.8 


0.08 


4.4 


0.20 


UM 139 


-1.6±1.9 


16 


6.2 


0.5 


0.03 


0.3 


0.02 


UM 141 


-0.4 ±0.1 


4 


0.7 


0.5 


0.12 


0.9 


0.25 


UM253 


-7.3 ±1.9 


114 


29.2 


0.7 


0.01 


2.0 


0.02 



M0~'* ergcm-2 s~' 

''Expected flux / minimum detectable flux 
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Fig. 2. — Same as Figure 1 for all objects observed with NICMOS camera 2. Each image is 3 x 3 arcsec. 
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Fig. 3.— 



Same as Figure 1 and Figure 2 for all objects observed with NICMOS camera 3. Each image is 6 x 6 arcsec. 



